Holographic QCD in medium: a bottom up approach 
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Abstract 

A holographic dual of hadrons at finite density is considered. We use the zero black- 
hole mass limit of Reisner-Nordstrom (RN) AdS background with hard wall to describe a 
confining background with finite quark density. We calculate density-dependence of meson 
masses and decay constants. In our model, pion decay constant and its velocity go down 
but all the meson masses go up as density grows. 
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1 Introduction 



Recently, motivated by the AdS/ CFT correspondence U, it has been shown that hadrons can 
be studied by constructing a 5 dimensional theory where chiral symmetry is imposed and the 
gluon dynamics is encoded in gravitational warping. This approach is generically called holo- 
graphic QCD H2H3HH. Five dimensional classical theory gives results of four dimensional 
hadron physics, fitting well both meson and baryon spectra. The scale in this model is in- 
troduced as cut off radius of the IR region. For the finite temperature case, one has to solve 
the wave functions in the black hole background and need to impose the infalling boundary 
conditions for the classical field at the horizon. The result shows that mesons are dissociated 
resulting in deconfinement at some high temperature 0. Confinement is described by the 
Hawking-Page transition where thermal AdS is the winning candidate in low temperature over 
the black hole background |6|. 

In this paper we are interested in nuclear matter where the baryonic density is high. One 
of the characteristic properties of nuclear matter is that it comes with others. No quarks are 
coming alone and nucleons in a single nucleus are coming in multitudes that are in interac- 
tions. While properties of hadrons, both mesons and baryons, in medium are fairly accurately 
understood up to the nuclear matter density thanks to both large amount of experiments per- 
formed since many decades and highly sophisticated many-body techniques, there is very little 
understanding of what happens when nuclear matter is squeezed beyond the density of nor- 
mal nuclear matter. This is principally because due to basic difficulty in dealing with fermions 
in the nonperturbative regime of QCD, there are no model-independent inputs available from 
theory and due to absence - up to date - of facilities probing high density, there is little guid- 
ance from experiments. 

The holographic approach to encoding the (baryon) chemical potential was suggested 
in and has been discussed by many authors HI . While the prescription to the geometry for 
the finite temperature has been very clear from the early days of AdS / CFT, it was far from the 
case for the corresponding geometry for the baryon charge. This is partly due to the fact that 
the very concept of the fundamental representation was introduced only by means of "probe 
branes." 

In Ref.0, it was suggested that if one introduces the bulk-filling flavor branes, the grav- 
ity back reaction of the baryonic matter can be easily encoded, allowing the charge of the RN 
AdS black hole to be identified with the fermionic charge as briefly explained in the next sec- 
tion. This back reaction should correspond to the attractive interaction between the gluon and 
fermion. The hydrodynamics and transport coefficients of this system have been analyzed in 
llTOll . However, what the low-temperature pair for the RN black hole which is the analogue of 
the thermal AdS as the pair of AdS black hole could be, has not been clear until very recently. 
In lfTT| , it was suggested that zero energy gravity solution is the resolution of the puzzle. There 
is a naked singularity, but it is hidden inside the hard wall. This background has the potential 
that can describe - for the first time - density dependence of a variety of observables in the 
confining phase. Especially, the phase diagram for the deconfinement phase transition and the 
p-meson mass were evaluated. At the present paper, we will make a quantitative study on 
other physical observables in finite baryon density, as well as correcting the numerical error 
for the p-meson mass in the section 4 of fTT| . Here, the masses and the decay constants of the 
vector mesons p, p' and the axial mesons it, a\ , a2 are calculated. We find that the meson masses 
typically increase at increasing density. Pion decay constant and its velocity go down but all 
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the meson mass goes up and as density grows. 



2 The Model 

To describe baryon density in a holographic setup, we will consider the RN AdS black hole 
geometry as background. To interpret the charge of the black hole as quark/baryon charge 
rather than the R-charge, we make use of a simple holographic model discussed by one of us [9] 
in which fermions are included as an local U(l) charge source on the bulk filling branes and 
the fermion back reaction on the metric is included by solving minimally coupled U(l) source 
to the gravity. The holographic model consists of N c D3 branes with Nf D7 branes touching 
the D3 branes and filling all of the AdSs space. The filling is optimal for D9 branes, but the 
differences between D7 and D9 are minimal in practical sense. The induced metric on D7 with 
zero current quark mass is identical to that of AdSs and the corresponding boundary gauge 
theory is SYMi + 3 with Nf massless flavors. We assume that the bulk filling brane tension does 
not modify the background metric due to its homogeneity while the fermion charge, which 
is dual to the local U(l) charge concentrated inside the black hole, requires the back-reaction 
of the metric. If we truncate the compact S 5 part and consider the gravity coupled with a 
local charge, the result is the well-known Reissner-Nordstrom AdS (RN-AdS) black hole metric 
coupled to a U(l) flavor bulk field sourced by the fermion charge on the black hole. The point 
of this construction is that it allows to interpret the local U(l) charge as the U(l) brane charge 
dual to the global fermion number rather than the R-charge. 

To describe the confining phase of this model we introduce a hard wall and set the mass 
of the black hole to zero, which corresponds to the zero temperature with finite density. The 
metric and electric potential of charged AdS black hole background are 

ds 2 = l l{-f{ z )dt 2 + dx 2 + ^ (1) 

f(z) = l + q 2 z\ A t = fi-Qz 2 

where I is the AdS radius. The black hole charge parameter q and the gauge charge Q will be 
related by the equation of motion it will be related to the quark number density. We put the 
fields in the region < z < z m and a hard wall at z = z rn . We can interpret this hard wall as 
a probe D7 brane in the background N c D3 brane [9J. This metric has no horizon hence has a 
naked singularity, which however is hidden inside the hard wall, which is typical in most of 
the confining background. We interpret this background as the low temperature pair of the RN 
black hole analogous to the thermal AdS as pair of the AdS black hole. Therefore we propose 
that it describes the confining phase in finite baryon density following MTT1IT21 . 
We start with the action (3) 



S = J d A xdzy/^gTr 



-\D M X\ 2 - M 2 X \X\ 2 - -^(L MN L MN + R M nR MN ) 



(2) 



We assume the metric is fixed and given by |(T]). X is the dual to the quark bilinear operator 
< qq > and M| = -3/1 2 according to the AdS/CFT dirctionary M\l 2 = A X (A X - 4). And the 
Tr is trace over the flavor group. The covariant derivative is 

D M X = d M X + iL M X -iXR M . (3) 
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If we write X = Se* 71 " 1 "" with real scalar S and adjoint ir a , the chiral symmetry breaking is ob- 
tained with nonzero vev of S and 7r a becomes the Nambu-Goldstone boson of chiral symmetry 
breaking. Here, M, N run over 0,1,2,3, z and fi, v run over 0,1,2,3. From L(eft) and R(ight) gauge 
fields, we form the vector and axial vector gauge fields 



V M 



-(Lm + Rm), A 



1 



M 



(Lm - Rm)- 



(4) 



Note that the length dimensions of the fields are [X] = L~ 3 / 2 , [V] = [A] = L" 1 , [tt] = [ip] = 
L°. Thanks to the rotational invariance of the boundary theory, i.e., SO(2), we can choose the 
direction of the wave propagating to be x 3 , so k^ = (uj, 0, 0, k). In the axial gauge, V z = A z = 0, 
the equations of motion of the vector fields are 








zf'(z) - f(z) 



1 



*/(*) 
k 2 



d z + 



u) 



2 1 



Vi(z) 



-d z + 
z f(z) 



V (z). 



(5) 



There is an extra dimension denoted by z, which gives rise to the towers of vector and axial 
vector mesons in 4 dimensions. Next we expand the 5 dimensional vector fields in orthonormal 
complete bases h n 's and a n 's 



V l (x,z)=^vt\*K(z), A i {x,z) = Y J Af{x)hi{z), A {x,z)=Y,A^ ] (x)ai(z). (6) 



We will consider only the time-like mass, so Vq can be ignored. The wavefunction should 
be normalized so as to give the proper four-dimensional kinetic term of the boundary gauge 
theory: 



I 



Zm dz {hX{z)y 



zf(z) 



1, I 



Zm dz (K(z)f 



zf(z) 



h I 



1. 



(7) 



In order to solve these second-order differential equations, we need two boundary conditions: 
we take the Dirichlet condition at z = 0, and Neunmann condition at z m , i.e., 



V(0) = 0, and V'(z n 
The equations for the axial-vector fields are 



0. 










zf'(z) - f(z) 



d z + 



zf(z) ~* ' f(zf 
A {z) 



2 v(zfl 
9t 



2;2l 



z 2 f(z) 



Ai{z) 



-d z — gl , 
z z 2 f{z) 



,v(z) 2 l 



2/2 



[d z ~ -d z ]<p - 9i^rj^{K + <p) 
2o 2 l 2 v{z) 2 f(z) a 



(8) 



(9) 
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where v(z) is the expectation value of the scalar S, (S) = \v(z), i=l,2 and the time-like mass is 
defined as to 2 = m 2 with k = 0. Note that A x 3 = d%ip is the longitudinal part of the axial field. 



Let us define g 2 = g 2 l, then all terms which contain g 2 v 2 l 2 are expressed as 

glv{z) 2 l 2 = g 2 v(z) 2 (10) 
where v(z) = l 3 ^ 2 v(z). The decay constant of the axial mesons and the pion are defined by 
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F 2 = — 

an n 2 



hf(z ) 



{F t,sy, 



(11) 
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3 Density Dependence of Physical Quantities 
3.1 Fixing the parameters 

To start with we fix the parameters of the model that can be obtained from what is available. 
The quark number density is N c times the baryon density i.e., n q = 3n# for iV c =3 and the 
nuclear saturation density is 0.16 fm~ 3 , so 

3 



B 



n 



47rr 3 
3n^ = 



0.16fm" 



1.28 x lO^GeV 3 



3.84 x 10- J GeV d 



(12) 



Normal nuclear matter is in the phase where chiral symmetry is spontaneously broken and 
color is confined. Hence the first-order phase transition from hadronic to quark-gluon phase 
predicted in the holographic models HJ [TTH must have a critical density above the normal nu- 
clear density. Quark number density n q is related to Q and q through k and 55 0: 



Q 



3 flf? 



2 / K 2 

We adopt the vacuum parameters fixed by p): 



n„ 



2 4q 

9s 



V& 2 

^5 



(13) 



Input 


value 


Input 


dill 


12tt 2 /N c 


Matching OPE of pQCD and hQCD 


m q 


2.0MeV 


m n and f n 


a 


(323MeV)' d 


m n and f w 




l/(323MeV) 


m p 



(14) 



remained is the gravitational constant n 2 . We determine k by comparing gravity with QCD 
thermodynamics: For QCD, the free energy is 



71 

45 



N 2 T\ 



while its gravity counter part in zero quark density is 

i 3 7r 4 r 4 

ads — q 7T ' 



(15) 



(16) 




By comparing these two, we get the relationship between k and N c 



2Ni 



2 ' 



For N c = 3, Nf =2 



Q 



2 12iV c 2 



457T 



iV„ / iV, 



11,: 



—$\ -^q ~ 0.0784831 g. 



(17) 



(18) 
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In this parametrization, one normal nuclear density is related to q ~ 0.0489278. Since different 
authors have different normalization conventions, here we give the summarized below. 



Coeff. 


PR[4,[13] 


BFB[9j 


Ours 


1 3 /k 2 


4N 2 /tt 2 


N 2 /4n 2 


2N 2 /45ir 2 




N c /12ir 2 


N c N f /4ir 2 


N c /12ir 2 




48iV c 


N e /N f 


8N C /15 


origin 


VVOPE 


SYM free energy 


QCD free energy 



(19) 



Note that in (4), Nf=2, the Nf factor is hidden in Trace operation. 
3.2 Chiral scalar 

The equation of motion (EOM) for the vacuum expectation value of the scalar field v(z) that 
comes from (X(z)) = v(z) = l 3 ^ 2 v(z) is 



zf'(z) - 3f(z) 



d z + 



z 2 f(z) 



v(z) = 0. 



The exact solution for v(z) is 



[\W ~q 2 z Q \ +az 3 2 F l Q, |, i, -q 2 z 6 



(20) 



(21) 



According to the AdS/ CFT dictionary, m q and a correspond to the current quark mass and 
chiral condensate (qq), respectively. Following |3]|, we take m q =2.29 MeV and a = (327MeV) 3 
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Figure 2: Left: First three axial-vector meson masses, Right: First three axial-vector meson 
decay constants \jF an in the range < q < q c ~ 2.5 no 



to reproduce the known values of pion mass and pion decay constant, and z m = (323MeV)~ 1 
to give the correct p meson mass 770 MeV at zero density. One of the disadvantage of hQCD 
is that there is no way to determine the chiral condensate. It should be provided as an input. 
In case of top down approach it can be easily calculated by studying the the embedding of the 
probe brane. 

3.3 Vector and Axial-vector mesons 

The spectrum of the vector mesons is obtained from the EOM (5). The procedure is the same 
as in I3j apart from the q dependence of the mass. For later discussions in confronting nature 
and comparing with the results of gauge theory models, it is important to understand that 
this is a "mean field" approximation in many-body theory language. In fig. deleft), the lowest 
p meson and its first two excitations are shown as a function of density. All the masses of 
vector mesons go up. The decay constants are plotted in figure [fright): they also increase at 
increasing density. 

The masses of axial vector mesons are also increased when density is increased. But 
unlike vector mesons, increasing mass of axial vector is almost ignorable for lowest mode, a\. 
See figure 2. 

3.4 Pion 

Our model shows that the space and time components of the pion decay constant, /* , /* , drop 
at increasing density. At the normal nuclear matter density, they are reduced to /* ~ 84.3 
MeV and /* ~ 91.2 MeV. The pion velocity defined as the ratio /£/ 'f* is shown in fig [4] Both 
fl and /* decrease as density increases, as expected from the chiral symmetry point of view. 
But the pion mass is increasing. One does not really know what happens to the pion mass in 
QCD with proper quark-mass terms. Experiments in deeply bound pionic atom indicate that 
the pion mass in medium actually goes up a few MeV at nuclear matter density. But since this 
increase persists in the chiral limit, this may be considered as a drawback of this model. Fixing 
this difficulty may require a nontrivial modification of the model as there seems to be no easy 
way to incorporate the "running" chiral condensate a. In the present treatment, a is introduced 
as an integration constant of the chiral scalar field v(z) ignorant of density of the system. 




III) 



Figure 3: Left: Pion decay constants versus quark number density. The dotted line is for f* and 
the solid line for /*., 




Figure 4: Left: Pion velocity Right: Pion mass 
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4 Thermodynamics 



It was shown in (11 J that a first-order phase transition between thermal charged AdS and RN 
AdS black hole takes place. The phase diagram without the calibrating the scale was drawn in 
the (fi, T), (Q, T) plane, but for different sets of parameters g\ , k 2 . In this section, we will briefly 
recast the work of lfTT| and draw the phase diagram with determined parameters in this work, 
which enables us to compare with real QCD phase diagram. 



4.1 Fixed chemical potential 

The difference of regularized action is given for the fixed chemical potential 



Srn 

Stc 

AS 



V 3 l 3 1(1 1 
k 2 T RN 
V 3 l* 1 



'-tc 



2n 2 [i 2 
3k 2 /x 2 



Srn — Si 



z: 



2g 2 l 2 %Tri 



tc 



k 2 T RN 



1 



3k 2 h 2 

rn 



24 2g 2 l 2 z 2 m 3g 2 l 2 z\ 



(22) 



the last equality comes from identification the Euclidean time periodicity at UV cutoff (e — > 0). 
This action is evaluated on physical state with Dirichlet boundary condition to the gauge field 
at UV, A T (e) = ifi. When this action difference is zero, there is first order phase transition 
which is identified to the de/ confinement transition in gauge theory side. Suppose that there 
is a critical point at z + = z c , then AS is zero at z + = z c . For the convenience, introduce 
dimensionless quantities, 



TrZf, 



then the action difference is rewritten 

V 3 l 3 1 



AS 



K 2 TRNzfn 



1 



2i 



rl 



+ 



k 2 a 2 

g 2 l 2 6z 2 + 



'Ml 



(23) 



(24) 



The temperature in dual hadronic background is identified with of RN AdS, 

1 



Trn 



■KZa 



1 



1 2 6 

2 q Z + 



and this q is identified to the chemical potential fi 

M =Q Z + = ^ q Z + 

The critical chemical potential and temperature in fixed chemical potential is 



(25) 



gji 2 3(1-2^) 

« 2 z 2 (9z 2 -2y 



1 

TTZ r 



~p? c 2k 2 „ 2 
^3g\ Zc 



(26) 



8 



T 








Figure 5: Phase diagram : number density (in no unit) vs T (in GeV unit). 



4.2 Fixed charge 



For the fixed charge case, we should add a boundary term to fix Q. This corresponds to impose 
Neunmann boundary condition to the at UV, then regularized action of RN AdS is 



qD 
D RN 



S» N + S b = S% N + \ [ d 4 xVGWn M A N F MN 
9r> JdM 



{0,0,0,0, 7 VTRK G^ = - A JJ(z) 



(27) 



where is the determinant of the boundary metric. This is [iN in statistical sense, so the 
grand potential is transformed to the Helmholtz free energy by adding this boundary action or 
imposing Neunmann boundary condition at UV cutoff. The action difference between RN AdS 
and thermal charged AdS is 



RN 



D tc 



AS 



V^J_ (I 

k 2 T RN 
V 3 l 3 1 

Srn — 



+ 



e 4 z 4 + 



Ak 2 Q 2 

w 

2k 2 Q 2 



1 



tc 



K 2 T R N 



2z% 



+ 



k 2 Q 2 
3g 2 P 



5zl - 2zi 



(28) 



the critical density and temperature is obtained after rescaling ((23]) 



Q c z 




2~z* 



z 2 l 



2) 

-2~zi 



2 5z? 



(29) 
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Here we set the IR scale z m = -^CeV 1 , then T c = 0.122 GeV. From the relation eq. flE), the 
critical density is 

n c = ^^q c ~ 1.1 x 10~ 2 GeV 3 ~ 2.86n (30) 
95 « 

The skyrmion approach tells that the critical density is almost two times of normal nuclear 
density HH). 



5 Discussion 



Using the zero black-hole mass limit of RN AdS background with its finite charge identified 
as the quark (or baryon) density, we studied a holographic dual of dense baryonic matter. The 
masses of the vector and axial-vector meson masses and their decay constants are obtained as 
a function of density. It is found that all the meson mass go up at increasing density. 

We have found that the pion decay constants decrease as density increases and the pion 
velocity given by the ratio of the space to time component of the pion decay constant also 
decreases. Chiral Lagrangian with pions alone would predict that the pion velocity goes to 
zero at T c . Our model has a glaring defect which needs to be rectified. The pion mass is found 
to increase more rapidly than the inverse of pion decay constant squared, thereby upseting 
in medium the GMOR relation valid in the vacuum. Also we can not calculate the density 
dependence of the chiral condensation from this model, unfortunately. The cause for these ill 
behaviors are not known at the moment. But we conjecture that it is caused by several elements 
that are left out in the model. Even the back reaction of the scalar fields which represent the 
effects of the gluon and chiral condensates has been ignored. We will report these issues in 
later works. 
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